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Cooperative phenomena arising due to the coupling of individual atoms via the radiation field 
are a cornerstone of modern quantum and optical physics. Recent experiments on x-ray quantum 
optics added a new twist to this line of research by exploiting superradiance in order to construct 
artificial quantum systems. However, so far, systematic approaches to deliberately design super¬ 
radiance properties are lacking, impeding the desired implementation of more advanced quantum 
optical schemes. Here, we develop an analytical framework for the engineering of single-photon su¬ 
perradiance in extended media applicable across the entire electromagnetic spectrum, and show how 
it can be used to tailor the properties of an artificial quantum system. This “reverse engineering” of 
superradiance not only provides an avenue towards non-linear and quantum mechanical phenomena 
at x-ray energies, but also leads to a unified view on and a better understanding of superradiance 
across different physical systems. 


A single atom coupled to an environment is usually subject to spontaneous emission and experiences a frequency 
shift referred to as the Lamb shift. In an aggregation of atoms coupled via the radiation field, collective effects can 
significantly alter the properties compared to a single emitter. For instance, this was realised by Dicke mm, who 
showed that N identical atoms confined to a volume much smaller than a wavelength cubed collectively behave as 
one “super atom”. This leads to exaggerated properties such as an acceleration of spontaneous decay by a factor of 
XDicke = Af (known as superradiance) or an enhanced frequency shift (sometimes also termed “collective Lamb shift” ) . 
Recently, also the correlated emission from extended ensembles of emitters has become the focus of experimental and 
theoretical mm investigations, where either the system size and/or the minimal interatomic distance a exceeds the 
scale of the characteristic wavelength Aq. The systems considered cover a wide range of possible realisations, including 
atoms near a nanofiber [5], thin vapor layers [7], cold atomic ensembles [SHU], or thin-film cavities with embedded 
Mossbauer nuclei in the realm of x-ray quantum optics [151 - 120] . 

The present work is motivated by the observation that in particular the latter experiments in the field of nuclear 
quantum optics exploited a deliberate control of superradiance properties, going beyond a mere characterisation. For 
instance, the observation of electromagnetically induced transparency at x-ray frequencies |13j was enabled by the 
engineering of two distinct ensembles with different superradiance properties in a single sample. Another example is the 
implementation of spontaneously generated coherences |14| . which relied on the realization of a spatially anisotropic 
electromagnetic environment via superradiance. In both cases, superradiance was employed to design an artificial 
quantum system, which in turn enabled the observation of the desired effect. 

This raises the question whether a systematic and constructive approach could be established to exploit superra¬ 
diance for the design of artificial quantum systems. Such design capabilities could overcome the limited resources 
accessible in state-of-the-art experiments, and thereby enable more advanced level schemes required, e.g., for the 
exploration of non-linear and quantum effects at x-ray energies. 

Here, we address this question by developing an analytical framework for superradiance in extended media en¬ 
compassing different system dimensionalities, interatomic couplings, and environments. As our main result, we then 
derive expressions describing how collective decay rates and frequency shifts can be controlled in extended media, and 
show how they can be used for the design of an artificial optical transition. 

We start with a single two-level atom (bare transition frequency uiq = cko, ko = 27r/Ao, c is the speed of light) which 
is embedded in an electromagnetic environment (e.g., free space) and is characterised by its spontaneous decay rate 
Re(Vo) = 7o (assuming Markovian reservoirs [21]). The coupling to the environment also results in a frequency shift 
Im(Vo)/2 = SujQ (single-atom Lamb shift). In the presence of an identical, second atom, photons can be exchanged 
between the two atoms. Due to irreversible loss to the reservoir, the inter-atomic coupling V^rj = Trirj + SiJwrjrj 
is complex [SHli]- Here, qrjrj ((^WriTj) represents the real-valued cross-damping (cross-coupling) term for two atoms 
located at positions and r^, respectively. Considering all pair-wise couplings in an ensemble of A^ ^ 1 atoms, we 
find 
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Figure 1. Design of an artificial optical transition through tailored superradiance. A d-dimensional lattice of atoms 
is embedded into an electromagnetic reservoir that mediates an inter-atomic coupling Vr oc l/r“, where atoms are separated 
by a distance r and the coefficient a characterises the distance-dependence (see eq. §). We show that the resulting collective 
eigenstates can be utilised for the implementation of an artificial transition with tunable decay rate and transition frequency. 


where E denotes the complex eigenenergy of the collective single-excitation atomic state I'F) = ‘PriCr^lO) 

and |0) signify the atomic raising operator for atom i and the vacuum state, respectively). Equation Q is valid for 
all dimensions d of the atomic arrangement and for all (physically reasonable) couplings 14.r - Collective decay rates 
and frequency shifts are obtained via E = —2Im(i?) and A = Re{E) — ojq, respectively 1^ . 

In Dicke’s small-volume limit, all atoms couple to each other with equal strength, leading to a collective decay rate 
r = Njq = XDicke7o ^nd a frequency shift A = XDicke<^wo with an enhancement factor Xpicke (see methods). To 
describe an extended sample, we consider ordered atomic arrangements, and focus on chains {d = 1), square lattices 
{d = 2), and simple cubic lattices {d = 3), see fig. [I The smallest inter-atomic distance is given by the lattice 
constant a. Such ordered arrays are naturally provided by crystalline samples (e. g., solid state targets employed in 
x-ray quantum optics HMH], optical lattices of atoms |26j . or atom-cavity networks m)- Furthermore, we consider 
a generic class of inter-atomic couplings 


K 

7o 


Ad sin^ 6 


^eikor 


(a > 0, e = ±), 


( 2 ) 


which depend on the distance r between atom pairs. Here, the coefficient a classifies the distance-dependence and 
Ad is a dimensionless coupling strength. We also assume the atomic dipole moments to be uniformly aligned along 
the Xs axis. This orientation dependence is taken into account by the angle 9 (see methods for further details). 
Since multiple terms of type ([^ can be accounted for by a linear combination, in particular also the three common 
implementations of three-dimensional free space [22] , atoms confined to two spatial dimensions |28j , or atoms coupled 
to a one-dimensional waveguide PU] are covered. The coupling parameters for these three examples are specified in 
table Note that in principle a can be artificially engineered and controlled as has recently been demonstrated at 
optical frequencies [10]. 


RESULTS AND DISCUSSION 

The solution of eigenproblem 0 (see methods) reveals that those eigenstates |k) whose wavevector’s magnitude 
matches the wavenumber set by the single atom transition, i. e., fc = |k| = ko, exhibit the maximum possible decay 
rate Tmax = Xmax7o if the constraint 

0 < a < ^ (3) 

is fulfilled. This criterion is a necessary condition for the emergence of superradiance and represents bounds on the 
allowed power laws of the coupling terms (exponent a in eq. 0 ) as a function of the lattice dimension d. For the 
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Table 1. Dimension-dependent quantities. The table summarises the quantities appearing in eqs. l|4|, ra, ig, and @- 
(191 as function of the system dimension d for the three considered example cases. Here, = arccos(kej; 3 /A:) denotes the angle 
between the eigenstate’s wavevector k and the *3 axis. 
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Figure 2. Collective decay rates and frequency shifts. Decay rates (black dashed curve) and frequency shifts (blue 
solid curve) as function of the wavenumber for a = (d — l)/2. The figure is valid independent of dimensionality and coupling 
type, due to the scaling of decay rate Xk = (Xfc — l)/(Xmax — 1), shift Aj, = [(Afe — 5a;o)/7o]/[e(Xmax — 1)] and wavenumber 
^ = (fc — ko)ahd\/N■ Note the offset h between the extrema of the frequency shift and the decay rate maximum. 


remainder, we assume that eq. ([^ is satisfied. The enhancement factor Xmax is (see table for quantities bd and c^, 
and methods for the prefactor Ld{a)) 
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In contrast to the maximum collective decay rate, we find that the collective frequency shift at fc = fco is always 
zero independent of the actual physical realisation. We thus conclude that the case of maximum superradiance is 
unsuitable for a control of both collective decay rates and frequency shifts. 

To circumvent this problem, we also consider states with wavenumbers around k = k^. Indeed, for a large but 
finite system, also states with a wavenumber close to /cq can exhibit an enhanced decay rate. We illustrate this for 
the most relevant case a = (d — l)/2 (which includes the three common implementations mentioned below eq. <§)■ 
For \k — fco|a ^ 1, we find 


— = Xfc - 1 + (Xmax - 1) • sinc(^), 
7o 

Afc - dwo Xmax - 1 C0S(5) - I 

70 2 ■ e 


( 6 ) 
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where ^ = {k — kQ)abd^/N and sinc(^) = sin(^)/^. Results are shown in fig.[^ scaled in such a way that they encompass 
different dimensions and coupling types. As mentioned before, those states which are maximally superradiant at 
k = ko do not exhibit a collective frequency shift. Rather, the frequency shift’s first two extrema around ko occur at 
wavenumbers k± = k^dih/{abd'^/N) (where h ~ 2.3311). This finding represents a unique feature as it is independent 
of the actual realisation and provides a signature suitable for a direct experimental test. 

Equations ([6 ) and ([^ also offer means to design an artificial optical transition with desired decay rate and frequency 
shift. In fact, the enhancement factor Xmax represents a characteristic scale for both decay rates and frequency shifts. 
As expected, we find that the particle number N and/or the sample volume V can be used to control Xmax- But 
additionally, eqs. <0 and 0 explain how the dimensionality d, the type of the inter-atomic coupling as described 
by a, as well as the coupling strength to the environment can be used to manipulate the enhancement factor. This is 
of particular relevance, since these parameters could also be tuned in situ |10[I30) . However, as mentioned previously, 
these quantities are not sufficient to change the ratio between decay rate and frequency shift. This only becomes 
possible by also controlling the wave number k (see fig. [^. Experimentally, the wavenumber could be adjusted via 
the excitation angle of the probing light field. 
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From a broader perspective, our results also enable us to understand how superradiant states from different reali¬ 
sations can be compared and categorised. This is important, e. g., if superradiant ensembles realised using different 
individual constituents are to be combined to an effective artificial quantum system. To this end, suppose that we 
can control the atom number and the volume such that N ^ N = JnN and V —V = /yV, respectively, where Jn 
and /v are arbitrary positive real numbers. Under this transformation, the enhancement factor changes as 


V ~ _ f dr^!U-<i-2Q) 

Xmax ^ Xmax — JN y JV Xmax • 


( 8 ) 


This behaviour allows us to classify superradiant states from systems with different dimensionality and types of 
coupling. For instance, we may say that two extended samples characterised by {d,a) and respectively, are 

similar if they satisfy the same transformation rule Q (leading to (a—l/2)/d = (o' —l/2)/(i'). As an example, if a one¬ 
dimensional system {d' = 1) should “imitate” the superradiant state from three-dimensional free space (d = 3, a = 1), 
the electromagnetic environment would have to be “engineered” m such that a' = 2/3. Similarly, if an extended 
sample should realise small-volume superradiance, transformation ([8) must reproduce the transformation of a Dicke 
system, which is simply XDicke XDicke with yoicke = N and yoicke = N = fNXuicke- Hence, a = (1 - d)/2 > 0, 
which reveals that only extended samples in one dimension (d = 1, otherwise we would have a < 0) can behave 


“Dicke-like”. 


In conclusion, we have studied single-photon superradiance in extended media, and showed how superradiance can 
be engineered in such a way that an artificial optical transition with tunable decay rate and level shift is realised. This 
result provides the basic building block for a systematic approach towards engineering advanced artificial quantum 
systems via superradiance by design. A promising avenue for future studies is the extension of our work to coupled 
sub ensembles with the goal to design artificial multi-level atoms [T5] . 


METHODS 


For an extended lattice, the plane wave ansatz for eq. y yields the eigenstates’ decay rates Fk = 

-2Im(i?k) and frequency shifts Ak = Re(Ak) — ojq as 


bk = 7o + Re[Id(k)], 

Ak = 5 ujq + ilm[Id(k)], 


Id(k) - ^ 


ikr 


(9) 

( 10 ) 

( 11 ) 


Here, r = {xi, ..., XdY' denotes a d-dimensional lattice vector with components Xi = arii, i = 1,... ,d, rii = — -^/N (2 + 
1,..., -sfN /2, and ^/N is even. Likewise, k = (fci,..., kd)'^ is the wavevector of the collective atomic excitation. 
The sum runs over all combinations of {rii} except ni = • • • = = 0 and the couplings depend on the distance 

r = |r| = a^rx\ + • • • + between atoms. We assume the atomic dipole moments to be uniformly aligned along the 
x-i axis (e.g., by applying a weak magnetic field). Thus, for d = 1,2 the distance vector r (in the X 1 -X 2 plane) is 
perpendicular to the dipole moments, and for d = 3 we have to take into account the polar angle 6 = &t:ccos{x^/ r). 
Furthermore, we make use of the assumptions N ^ \ (many atoms) and k^a > 1 (extended sample). The decay rate 
eq. © can also be rewritten in terms of the enhancement factor 


_ Fk ^ , Re[Id(k)] 

Xk = — = 1 H- 

7o 7o 


( 12 ) 


To arrive at the final expressions eqs. 0-0, we further manipulate eqs. 0 -([I^ as follows. In this paper, we 
focus on the system’s eigenstates and—to keep the analysis general—do not consider geometric details or questions 
of how to excite and probe the system since such details vary from experiment to experiment. Around fc = |k| = k^, 
we can utilise a continuum formulation, rewrite the lattice sums in eq. (11) into an integral, and perform the angular 
integration for couplings of type 0 (see supplementary information for further technical details of the calculation). 
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The dimension-dependent quantities Ad^ bd, Cd, and gd are listed in table (for instance, Ad is real for d — 1,2 
and purely imaginary for d — 3). Note that the factor exp(±ifcor) from eq. (|^ in the eigenproblem Q can be 
understood as a radial translation in wavenumber space. In the shifted frame, a long-wavelength limit of the collective 
atomic excitation (which can be accounted for by a continuum description) corresponds to A: —>■ fep. This continuum 
formulation is applicable in the range |A: — fco| tt/ abd\fN. Further, in eq. (|^, we have not included exponential 
damping of the form exp(—fepr/^), where ^ denotes a dimensionless absorption length that, for instance, empirically 
accounts for material imperfections. Such a damping factor in the integral in eq. (14) would lead to a broadening and 


modification of the k = /cQ-criterion for maximal superradiance, going beyond the scope of this paper. Details on the 
calculation of the integrals in eq. (14) can be found in the supplementary material. 

The maximum enhancement factor (|^ can be cast into the equivalent forms (V = Na‘^ denotes the sample volume, 
p = N/V is the number density, and Xmax — 1 — Xmax since N ^ 1) 
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Which formulation to choose from eqs. (16)-(18) depends on which quantities can be controlled in an experiment. 


If for small volumes the length scale set by the inter-atomic distance a is effectively eliminated from the single¬ 
excitation eigenproblem 0 (possibly neglecting divergent contributions to the inter-atomic coupling [51 |3T]), all 
atoms couple to each other with equal strength Vq. The resulting equation E — ujq = —(i/2)Vo T’rj/T’r, (which 

must hold for all r^) yields a maximal decay rate for a spatially constant wavefunction with equal relative phase 
between all atom pairs, representing the maximally symmetric Dicke state. For this state, F = —2Im(i?) = Njq and 
A = Ke{E) — ujQ = NSujq- 
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SUPPLEMENTARY INFORMATION 

In this supplementary information, we provide technical details on how to rewrite the lattice sums into an integral 
and on how to ultimately perform the integration. 


I. TdiX) 


The quantity 
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where the sums run over all combinations of {rii} except ni = • • • = = 0, can be rewritten into an integral 


Xd(k) 


Ad^ „ 

^>'(.■,9)11 

i=i 


^-kjXj 


( 20 ) 


( 21 ) 


For d = 1,2, V{r,0) = fr, whereas for d = 3, V(r,9) = sin^ 0 fr {fr = Ad eTq){eikor)/{kor)°‘)■ Here, dAx signifies 
the d dimensional infinitesimal volume element and the integration is over all space except for a region with radius a 
around the origin (d^cc = dx, d^a; = rdrd(p, d^cc = A sin0drd(/?d0). 

Explicitly, for d = 1, 
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where N 2 is chosen such that the integration area covers N atoms, i.e., = N. For d = 3 
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X Jo ^sindy ^kl + 


(24) 


where (47r/3)(A3)^ = N, Jo(-) signifies the zeroth-order Bessel function of first kind, and sinc(x) = sin(x)/x. The 
integration over 9 can be done as follows. Upon definining Im = /J^ ddsin™ dexp)—ifcj_r cosd) Jo(A:||r sind), kx = k^, 
and fc|| = \/kl + we have the relation 
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To simplify the integration needed for Ii, we can choose a coordinate system in which either fc_L = 0 or fey = 0 (it can 
be shown that Ii does not depend on the orientation of k), yielding 
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= 2 sine 



(26) 


Finally, la = sin^ r? • 2sinc(fcr) + 0[{kr) and therefore 

dTT n 

X 3 = —^ sin^ r? / dr r^/i.sinc(fcr), (27) 

® Ja 

where r? denotes the angle between the eigenstate’s wavevector k and the 2 axis. Here, we have only taken into 
account the asymptotic leading order term (with respect to kr). Other terms can be accounted for by means of 
different coefficients a (see main text). Introducing the abbreviations used in the paper and performing a variable 
substituion, we finally arrive at the integrals Jd{k) (Eq. (14) in the main text). 


II. Jd(k) 


We now proceed with the radial integration 

koabd vW 

Jdik) = J dr] e^^’^gdikkg^r]) 77 ^. (28) 
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Note that for d = 2 the integration kernel is actually given by Jo{kkQ^ 7 ]). However, already at the lower integration 
limit, we can use the asymptotic form Jo{ka) « yj2/'KCOs{ka — TTl4)l'/ka since in an extended sample k^a > 1 
and only the wavenumbers k around fco are relevant. Furthermore, by an additional substituion of the integration 
variable, we shift the 7 r /4 shift to the argument to the exponential (which we can account for by means of appropriate 
prefactors), 77 + 7 r /4 « 77 , and the integration limits can also approximately remain unchanged. 

The possible combinations in the integrand we need to consider are 
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(31) 


(32) 


A. Jee 


Utilizing a computer algebra system, we find that 


Jni 


i/3+1 


(^) (sgn(r-)) 

\^o / 

X (^r(l + /3, -ir_a) 

- r(l +/3,-ii’-aA^d) 


(33) 
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- (-l)^r(l + /3, ir_a) 

+ (-l)^r(l+/3,ir_a7V')) 



X ^r(l + j3, — ir+a) 

- r(l +/3,-ir+a7Vrf) 

- (-l)'^r(l +/3,ir+o) 

+ (-l)^r(l + /3,ir+a7V^)) , 


where r± = fc ± fco and r(a, z) = dt ^e * signifies the incomplete Gamma function. The asymptotic form for 
1 and r_ 7 ^ 0 reads 


JcciJ^ 7^ ^ 0 ) 


jd+l 


' ^sgn(r_))- 2 ^ 

V ^0 / 

X (^r(l + /3, -ir_a) 

-(-l)^r(l + / 3 ,m-a) 

-2i(-l)/5(ir_aiV')''sin(r_aiV')) 



X (^r(l + /3,-ir+a) 

- (-l)^r(l + /3,m+a) 

- 2 i(-l)^(ir+aiV^)^sin(r+aA^rf)) . 


The dominant terms for N'^ ^ 1 in this expression are oc if /3 > 0. 

York^ko {N'^ = hd<fN->l), we arrive at 


Jcc{k —)■ fco) — 


2(1 + /?) 
{bdY^^ 
2 ( 1 + / 3 ) 


{k,af+^ {N'df^^ 


Here, the dominant terms for N'^ ^ 1 are (x if /3 > —1. 

For the case /3 = 0, the explicit expression for the integral reads 


Cl 


13=0 


koa 

^ L 


N'd {smc{r-aN'd) + smc{r+aN'd)) 
— sinc(r+a) — sinc(r_a)j 


|r_|a<l, Ar>l /coa , . , 

- ^^dSinc(r_alVrf) 

fc^o kpa , kpa 1 

-Y^d. = -Y^dN-i . 


(34) 


(35) 

(36) 


(37) 


In the second step, we focus on wavenumbers k around kp (i. e., small |r_ |a, for which the terms with r+ are negligible). 
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B. Je 


Similarly, 


Jsc - 


i/3+2 


(sgn(r-)) 

\ko J 

X ^ - r(l + /3, -ir_o) 

+ r(l + ^,-ir.alVrf) 

- (-l)^r(l + /?,ir_o) 

+ (_l)/3r(l + ^,ir_aiV')) 

(r(l +/3,-ir+a) 


+ 


- 1-/3 


-r(l + /3,-ir+olVrf) 

+ (-l)^r(l + /?,ir+a) 

-(-l)^r(l + /?,ir+alV^)) 


The asymptotic form for ^ 1 but r_ 7 ^ 0 reads 

i/3+2 

Jsc{k 7 ^ fco) — i ■ 


(38) 


X ^ - r(l + /3, -ir_a) 

- (-l)^r(l + /?,ir_a) 

+ 2(—l)^(ir_a-/V^)^ cos{r-aN'^) 

- 1-/3 

(r(l +/3,-ir+a) 

_(_l)/3r(l + /?,ir+a) 

- 2(-l)^(ir+aA^d)^ cos(r+alVrf) 

In contrast to the integral Jcc, the limit /c —> fco does not yield a scaling oc The case /3 = 0 reads 


(39) 


3=0 fcoo 


cos{r-aN^) cos{r+aN'^) 


T-a r+a 

cos(r_a) cos(r+a) 
r_a r+a 

|r_|a<i koa cos(r_aiV^) — 1 
2 r_a 


(40) 


C. J„ 


We can rewrite 


koabd vW 


dr-R - 


d ?7 cos{r]) sin(fcfcQ ^i/) 


fco a. 

fc 


/3+1 


fca/>d 


d ?7 sin(? 7 ) cos(fcofc ^ 17 ) 17 ^ 


(41) 
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This is in essence just the integral Jsc with k and ko interchanged. In particular, for (3 = 0 , Jcs = —Jsi 


D. Js. 


For k ^ ko, we can write 


koabd vdv 


Jss{k ko) = j dr] (1 - cos^ t]) 


2^) ^/3+l 


= (koabd^)^^^ -(A:oa)^+^ - Jcc(A: ^ fco) 


2J„c{k—^ko) 
— Jcc{k —)■ ko) ■ 

For (3 = 0 , the explicit expression reads 

/3=o koa r , 


N'd (sinc(r_alVrf) - sincir+aN'j)) 
+ sinc(r+a) — sinc(r_a) 


|r_|a<l. Ar>l fcoa , . , 

- ^^dSinc(r_aiVrf). 



